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FOREWORD 


This  memorandum  represents  a portion  of  work  being  done 
in  support  of  SF  101  03  16,  Task  8132  NUWC  Problem  Number 
Elll  by  the  Sonar  Signal  Recognition  Division. 

The  information  contained  is  considered  useful  to  those 
groups  engaged  in  analysis  and  simulation  using  wideband  signals. 

Appendix  A was  developed  by  L.  R.  Weill  and  was  previously 
presented  with  experimental  verification  in  "Analysis  of  Active- 
Sonar  Signals  (U),"  NUWC  TP  7,  Sept  1967,  CONFIDENTIAL  by 


D.  G.  Olson  and  L.R.  Weill. 

Only  limited  distribution  is  contemplated. 


I 

I 


i 


Q U 


INTRODUCTION 

Sn  many  circumstances  an  investigator  wishes  to  use  digital 
computers  for  simulating  sonar  receivers.  Quite  often  the  pro- 
posed implementation  of  beam-formers  or  matched-filter  re- 
ceivers will  differ  from  analytically-derived  filter  characteristics. 
In  order  to  evaluate  different  implementations  of  receivers  with 
respect  to  each  other  and  the  theoretieally  ideal  version,  an 
investigator  must  use  a valid  representation  of  a sonar  transmit 
signal  in  the  computer  simulation. 

Several  references  41=4f1iave  described  some  of  the  differ- 
ences between  narrowband  signal  representation  and  wideband 
signal  representation.  However,  much  of  the  signal  processing 
literature  is  directed  to  the  radar  problem  which  is  usually  a 
narrowband  signal  case.  Therefore,  this  paper  will  collect  the 
comments  of  several  authors  to  demonstrate  the  differences  of 
narrowband  and  wideband  assumptions  and  their  effects  on  the 
concepts  of  complex  signal  notation,  complex  modulation,  and 
signal  envelope . 

S, 

REAL  SIGNAL  RElPRESENTATlON 

‘\ 

Any  transmittable  radar  or  sonar  signal  belongs  to  a class 
of  functions  describable  as 


s(t)  ^ a(t)  cos  [2TTfQ  t + !p(t)] 


(1) 


Whether  this  equation  provides  a particularly  useful  representa- 
tion, however,  depends  on  the  problem.  In  equation  (l)a(t)is  an 
amplitude  modulation  function,  fy  is  a carrier  frequency  and ^(t)  is 
a time-varying  phase  modulation  function.  The  argument  of  the 
cosine  is  the  signal  phase  and  the  derivative  of  the  signal  phase  is 
defined  as  the  instantaneous  frequency 


" 'o  2tt  “dT 


(2) 


For  a given  real  signal  s(t),  the  choice  of  fy  is  arbitrary  since  a 
suitable  ^{l)  can  be  chosen  to  yield  tlie  given  value  of  s(th  There 
is,  however,  a physically  meaningful  waj'  of  defining  the  carrier 
frequency  of  a modulated  signal. 

Before  proceeding,  the  use  of  the  Fourier  transform  should 
be  introduced  where 


1 


1 


s(t)  = 

S(f)  clf 

J~oa 

H signal 

(3) 

S(f)  = 

f s(t)  e‘j2Ttft  fit 

J — oc 

= frequency 
spectrum 

(4) 

and 

r [s(t)]2  dt  = r 

oo  J 

|S(f)l^  df  = signal  energy 

(5) 

The  frequency  spectrum  given  by  the  Fourier  integral  is  a com- 
plex function  of  f,  and  extends  over  all  positive  and  negative  fre- 
quencies. Since  the  transmitted  signal  s(t)  is  real,  s(t)  = s * (t) 
and  S(-f)  = S * (f)  , where  the  asterisk  indicates  the  complex 
C'-.ijugate.  Now  the  carrier  frequency  can  be  defined  as  the 
first  moment  of  one-half  the  energy  spectrum 


(6) 


Since  SiOextends  over  positive  and  negative  frequencies,  it  is 
called  the  "double-sided"  spectrum.  The  spectrum  of  a real 
signal  is  even  in  amplitude  and  odd  in  phase. 


where  S(f)  = lS(f)l  (7) 

with  |S(f)|  = amplitude  of  the  frequency  spectrum  (8) 

and  = phase  of  the  frequency  spectrum 


Because  one-hrdf  the  energy  is  contained  in  the  positive  spectrum, 
the  carrier  frequency  was  defined  using  E '2  in  equation  (6).  A 
pictorial  representation  of  a narrowband  signal  is  given  in  figure  1. 
While  the  first  signal  amplitude  spectrum  in  figure  1 represents  a 
narrowband  signal  with  no  frequencies  from  the  right-half  of  the 
amplitude  spectrum  spilling  over  into  the  negative  frequency  por- 
tion, the  second  amplitude  spectrum  represents  the  wide  band  case. 
All  of  the  previous  equations  apply  for  either  narrowband  or  wide- 
band signals.  Now  the  following  approximations  will  show  the 
utility  of  equation  (1)  for  representing  narrowband  signals.  Using 
equation  (5)  with  (1)  gives 

E - j cos  I27rf^,  t I clt 

- H- (t' clt  » —y  a- (t)  cos  2 [2Tr  fy  t • 9 (t)]  ell  (10) 
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Amplitude  Spectra  vs  Frequency  For  o Wideband  Signol 
Figure  1. 


f 


In  the  case  of  a narrowband  signal,  both  a(t)  and  9 (t)  vary 
slowly  in  comparison  to  cos  2n  t.  The  second  integral  on  the 
right  hand  side  of  (10)  has  a product  term  formed  from  the  two 
functions  as  shown  in  figure  2 such  that  the  integral  of  the  product 
is  virtually  zero.  Thus,  as  a consequence  of  narrow  bandwidth 


1 /■« 


2 (t)  dt 


(11) 


This  approximation  greatly  simplifies  calculations  as  it  permits 
the  substitution  of  the  amplitude  modulation  function  for  the  signal. 

COMPLEX  SIGNAL  REPRESENTATION 

As  a further  simplification  for  theoretical  analysis  involving 
linear  operations,  the  real  signal  is  often  expressed  as  the  real 
part  of  a complex  waveform  (t) 


s(t)  Re  !<,'/(t)l 


(12) 


where  Re)  I means  the  real  part  of  ) I-  The  analyst  replaces  sit) 
with  ^ (t)  for  these  linear  operations  and  then  takes  the  real  part 
of  the  calculation  for  his  answer.  First,  rewrite  equation  (3)  as 


s(t)=P  S(f)  ei2"ft  df  - S(f)  cj2Trft  df 

•/  0 J —eo 

and  using  S(-f)  = S«(n  for  real  signals 

s(t)  = P S(f)  ci2TTft  df  . s * (f)  e-j2TTft  (If 

-'O  Jo 

which  becomes 

s(t)  = 2 S(f)  ei2Trf(  df] 


(13) 


(14) 


(15) 


Since  the  negative  frequencies  simply  mirror  the  positive  fre- 
quencies in  complex  conjugate  form  for  real  signals,  the  negative 
frequencies  can  be  omitted.  Equation  (16)  shows  that  s(t) can  be 
regarded  as  the  real  part  of  a complex  signal  whose  frequency 
spectrum  is  twice  that  of  tlie  real  signal  for  positive  frequencies 
and  zero  for  negative  frequencies.  Thcrefoi'e,  using  >1’ (f)as  the 
complex  signal  frequency  spectrum 


t (f)  - 2S(|) 
0 


f > 0 
f < 0 


(16) 
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cos  2 [2iTfQt+  9(t)] 


a^{t)  cos  2 [2TifQti  9(t)] 


Figure  2.  The  integral  of  the  product  a^(t)  cos  2 [2TTfQ  t i 9(t)]  for  a slowly  vary- 
ing o(t)  and  9(t)  will  nearly  equal  zero  as  the  contributions  above  and  below  the 
axis  will  practicolly  cancel  each  other. 


and  'I'  (f)can  be  called  a "single-sided”  spectrum. 

CONSTRUCTING  THE  ANALYTIC  SIGNAL-  or  COMPLEX 
WAVEFORM 


While  the  frequency  spectrum  of  the  desired  complex  signal, 
'I'  (f) , has  been  stated,  an  analyst  still  needs  to  find  a suitable  i>  (t) 
corresponding  to  a given  s(t).  Consider  now  the  characteristics 
of  the  Hilbert  transform  defined  as 


s(t)  = 


s(t) 

(t-T) 


dT 


(17) 


where  the  principal  value  is  intended,  and  s(tHs  a real-valued 
function  on  -~  < t < +=».  References  (5-6)  show  several  character- 
istics of  this  transform.  From  equation  (17)  the  Hilbert  transform 
is  seen  to  bo  a convolution  of  s(t)  with  the  function  -(ttI)''  , written 
as  s(t)  * (irt)'! . Remembering  the  property  of  the  Fourier  trai\s- 
form  where 


g(t)  * h(t)  G(f)  11(f) 


(18) 


or  convolution  in  the  time  domain  transforms  to  multiplication  in 
the  frequency  domain,  and  using  the  Fourier  transform  of  (nt)’' 
as-jsgnf,the  spectrum  of  s(t)  is 


-j  S(f) 

f > 0 

0 

f = 0 , 

(19) 

jS(f) 

f < 0 

If  we  choose  ^(r(t)as  the  pre-envelope  discussed  in  reference  (6) 
where 


tjj  (t)  = s(t)  + js(t) 


(20) 


the  Fourier  transform 

>1'  (f)  = S(f) 


-jS(f) 

f > 

0 

f = 

jS(f) 

f < 

2S(f) 

f > 

8(0 

f = 

0 

f V 

Thus,  by  finding  s(t)from  s(t)  and  writing 


(21) 
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ijj  (t)  = s(t)  + j s(t) 


(20) 


a complex  signal  can  be  formed  from  a real  signal  with  the  fre- 
quency domain  behavior  shown  in  equation  (16). 

Now  the  Hilbert  transform  of  cost  is  sin  t.  Reference  (5) 
lists  other  transform  pairs,  but  defines  the  Hilbert  transform  as 


X T"" 

' J-o^  T-t 


which  is  the  negative  of  equation  (17).  Therefore,  the  transforms 
and  formulas  in  reference  (5)  must  be  multiplied  by  (-1)  for  equiva- 
lence with  the  expressions  in  this  paper.  Thus,  considering  a 
narrowband  signal  with  a slowly  varying  phase  modulation  function 
9(1), 

Jjji  i cos  [2tt  fp  t + 9 (t)  ] i = sin  [ 2 tt  t + 9 (t)  ] (23) 

where  Jji.!  1 indicates  the  Hilbert  transformation.  For  this  narrow- 
band  case  the  complex  waveform  corresponding  to  s(t)  becomes  to  a 
close  approximation 

= ait)  cos  [2tt  fyt  f 9(t)]  + j a(t)  sin  [2  7Tfot  1 9(1)]  (24) 

Stating  equation  (24)  as  an  equality  and  simplifying  the  notation  by 
use  of  tile  exponential,  we  get 

, ,,,  j[2Trfot  + 

(25) 

= Hit) 

where  uit)  is  a complex  modulation  term 

Hit)  = ait)  c>‘P'*'  (26) 

consisting  of  the  real  amplitude  modulation  function  and  a complex 
phase  function.  The  evelope  of  the  signal  is 

|</it)i  = -t  s^it)  (27) 


which  gives  U'/it)!  = ait) 


for  this  narrowband  case.  Wliile  equations  (20)  and  (27)  are  per- 
fectly general  and  follow  reference  (6),  the  results  shown  in  equa- 
tions (25)  and  (28)  are  valid  approximations  only  for  narrowband 
signals. 
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While 


3jlj  tcos  (2TrfQt  ^ 9 1 = sin  (2TrfQt  4 9)  ^29) 

for  a particular  value  of  9, 

?jji  I cos  [2  nf^t  4 <p(t)  ] i sin  [ 2iTfyt  4 (p(t)  ] (30) 

in  general  and  the  approximation  of  equation  (23)  only  holds  in  the 
narrowband  case.  For  wideband  signals  the  Hilbert  transform 
could  be  obtained  by  expressing 


s(t)  = alt)  cos  [2711^1  4 (p(t)l  (1) 

as  a Fourier  series  and  then  transforming  each  term  by  the 
addition  theorem  [reference  (5)] 

?Hi  •Wt).4  g(t)!  = lh(t)!4[Fjj;  jg(t)i  (31) 

In  summary,  the  representation  of  any  sonar  or  radar  trans- 
mit waveform  by  the  real  signal 

s(t)  ^ alt)  cos  l2TTfyt  4 9(t)l 

is  completely  general.  In  the  case  of  narrowband  signals  - which 
implies  a slowly  varying  alt)  and  9lt)  - the  complex  signal  u It) can 
be  readily  written  as 


.,)«)-  ult)ei2’^<'ot 


(25) 


and  the  evelope  as 


10  It);  = alt) 


(28) 


which  is  tlic  amplitude  modulation  fimction.  In  the  wideband  case, 
finding  the  Hilbert  transform  sit)  mai'  be  more  difficult  than  simply 
using  the  real  signal  i-epresentation  of  equation  (1)  for  analysis 
purposes.  For  conijuiter  simulation  purposes  which  use  time- 
series  samples  of  the  sonar  signal,  there  is  no  need  to  consider 
the  complex  signal  form  anjavay.  Appendix  A will  outline  a direct 
computational  mctliod  of  reconstructing  tlie  envelope  defined  as 

Hit)  - >/s -1 1 ' ‘ s" It ) (32) 

from  time-series  samples  of  the  real  waveform  sit) 
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USING  COMPLEX  NOTATION  FOR  WIDEBAND  SIGNALS 


The  preceding  section  discussed  the  difficulties  of  fiiuUng  a 
complex  signal  corresponding  to  a particular  wideband  real  signal 
specified  in  the  time  domain.  If  an  anal3'st  wishes  to  study  the 
sonar  problem  in  the  frequency  domain,  however,  he  can  define 
the  transmit  signal  as  a modulated  carrier.  The  approach  is  dis- 
cussed in  references  (3)  and  (4). 

Using  the  Fourier  transform  relations  of  equations  (3)  and 
(4)  and  the  consequence  of  s(t)  being  real  as  shown  in  equation  (15), 
we  can  write 


s(t)=2RcjJ^’°  S(f)ej2^''Mf 

= 2Rejej2^fo‘  f)ei2’^f‘df 


(33) 


If  we  define  a function 

/oo 

S(fo  ^ f)  eJ2TTftcif 

"^0 


(34) 


we  can  then  express  the  real  signal  as 

s(t)=2  RolV(t)ej2TTfoti 


(35) 


which  consists  of  a complex  modulation  function  and  the  ct  rrier 
frequency  term  in  the  exponential.  The  carrier  frequency  term 
is  arbitrary,  but  the  fimction  \’(t)  becomes  unique  once  f^  is 
chosen.  Naturally,  can  be  chosen  as  shown  in  equation  (6). 

To  further  examine  the  properties  of  \'(t)  , consider  the 
absolute  value  after  a change  of  variable 


such  that 


so 


fo  + f=f 

(36) 

df=  df' 

(37) 

V(t)  : 

=2  r 

Jo 

1 

(38) 

= 2c 

i2’Tfot  S(f')c)2’^f'*cip 

Jo 

(36) 

V(t) 

-I2  f"  S(r')oi2^i'\lp 

1 Jo 

(40) 
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Now  consider  the  method  of  computing  the  envelope  using  the  Hilbert 
transform  where 

ld(t)i  = ls(t)  + iSt)l  (41) 

Using  the  addition  theorems  of  the  Fourier  and  Hilbert 
transformations 

ld(t);  = IS'-I  IS(f)l  * ri-M-isgnfS(f)li 

= S(f)  cj2Trftdf  , j -jsgllfS(f)  ej2Trftjc 

= [j’~  S(f)  0>2^f‘clf  H S(f)oi2’Tftdf 

4 r°  S(f)ej2^f‘df-  S(f)  ej2TTftdfl 
J — oo  J —oo  » 

I/*  oo 

2 / S(f)ei2’Tftdf 

Jo  (42) 

where 7"'^l~jsgnfS(f)l=- St),  since  -jsgiifS(f)is  the  Fourier  transform  of 
s(t)  as  discussed  previously.  Clearly,  equations  (40)  and  (42)  are 
equivalent.  Therefoi’e,  iVtt)!  represents  the  envelope  of  the  sigrial 
and 


lV(t)l  = ls(t)  + js(t)l 


(43) 


Thus,  if  an  investigator  wishes  to  e.xpress  a wideband  signal  as  a 
modulated  carrier,  he  must  first  select  a carrier  frequency 
and  then  compute 


/oo 

S(fo  t f)  ej24Tft 

which  is  related  to  the  real  sigmal  as 


df 


s(t)  - 2HolV(t) 


(34) 


(35) 


This  process  does  not  require  the  use  of  the  Hilbert  transform  of 
the  real  signal.  Only  the  investigator  with  a particular  problem 
can  select  the  most  advantageous  wideband  signal  I’cprcsentation 


s(t)  a(t)  cos  l2TTf' ,t  ‘ 9(t'l 


(1) 
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APPENDIX 

If  a real  analog  signal  is  properly  sampled,  the  signal 
waveform  can  be  reconstructed.  If  the  transmitted  signal  s(t)  is  a 
real  function'defined  on  — ^ < t < oo  and 

(1)  s(t)  has  a Fourier  transform  S(f)  such  that 
sup  t f lS(f)  / 01  = fjj 


is  finite 


(2)  the  sequence  of  sample  values 
— f(-nAt),— , f(0),  f(At),  — f(nAt) 


is  known  for  all  integers  u 


then  s(t)  can  be  determined  exactly  from  the  sample  values  alone 
at  all  t where  f is  continuous.  Specifically,  the  reconstructed 
signal  s(t)  = s(t) 


when 


where 


.=  2 s(r 

n=  - 9c 

'"(■A  -") 

(m'") 

n - 0 


s(nA  t) 


"it  -") 


(A-1) 


(A-2) 


(A- 3) 


This  theorem  is  proven  in  reference  (7).  This  treats  the  base  band 
signal,  or  a signal  where  the  highest  frequency  f|,  is  known  rather 
than  the  actual  signal  bandwidth  W = f„  - . The  choice  of  sampling 

interval  At  when  both  fj,  and  f,  are  luiown,  and  when  the  first  two 
conditions  leading  to  (A-l)  are  satisfied  in  addition  to  the  conditions 

'■>“  Mf'o 


(2)  the  sampling  intci’val  At  satisfies  one  of  the  following: 


A-l 


K-1  K 

2f,,  " " 21,, 

K-2  . K-1 

2f.  ‘ ^ 2f,. 


if  "-r 

'^‘l  'h 


0 <At  < 


for  K the  largest  integer  such  that  K 


thens(t)  = s(t)  where 


s(t)  = At 


^s(n 


sin  27Tf„  (t-nAt)  sin2TTf,  (t-iiAt) 
2uf„(t-nAt)  “^^l^'2TrfL(t~nAt) 


as  shown  in  reference  (8) . 

Now  consider  the  Hilbert  transform  of  s(t) . The  Hilbert 
transform  of 


cos  x-1 


by  equation  (22) 


1 - cos  X 


by  equation  (17) 


1 2siM2(^) 

1 - cos  X \ 2/ 


(A-5) 


we  can  write  the  Hilbert  transform  of  the  two  reconstruction  form- 
ulas as 


s(t)  ==  > s(nAt) 


(ii  ~")l 

”(r,  -") 


(A-6) 


A-2 


sin2[Trf||  (t-iiAt)] 
2TTf„(t-nAt) 

sin^  “ iiAt)] 

^ 2 Trfj^  (t  - iiAt) 

Using  a sampling  interval  (a)  such  that 

t = ga,  g 0,  1,  2,  N-1 

a simplified  expression  for  computing  the  envelope 
R(ga)  ==  |S(ga)| 

2 sin^  y (g-n) 

TT(g-ll) 


and  for  the  other  case  of  known  bandwidth 


